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Abstract. The paper is devoted to optimization of resonances in a 1-D open optical cavity. The 
cavity's structure is represented by its dielectric permittivity function e(s). It is assumed that e(s) 
takes values in the range 1 < S\ < e(s) < e 2 . The problem is to design, for a given (real) frequency 
a, a cavity having a resonance with the minimal possible decay rate. Restricting ourselves to 
resonances of a given frequency a, we define cavities and resonant modes with locally extremal 
decay rate, and then study their properties. We show that such locally extremal cavities are 1-D 
photonic crystals consisting of alternating layers of two materials with extreme allowed dielectric 
permittivities e\ and e 2 . To find thicknesses of these layers, a nonlinear eigenvalue problem for 
locally extremal resonant modes is derived. It occurs that coordinates of interface planes between 
the layers can be expressed via arg-function of corresponding modes. As a result, the question of 
minimization of the decay rate is reduced to a four-dimensional problem of finding the zeroes of a 
function of two variables. 

Key words: photonic crystal, high Q-factor resonator, quasi-normal eigenvalue optimization, non- 
linear eigenvalue 
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1 Introduction 

A leaky optical cavity (or resonator) is a region of space within which the electromagnetic field is 
well confined, but not completely confined. Because of the leakage, each standing electromagnetic 
wave in the cavity decays exponentially in time. The rate of energy loss and the frequency of these 
eigenoscillations can be characterized by the corresponding complex eigenvalue u, which is called 
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a quasi-normal (QN) eigenvalue or a resonance (the latter name is used in Quantum Mechanics 
and, sometimes, Engineering). The imaginary part Imw of the QN eigenvalue corresponds to the 
decay rate of the standing wave, the real part Re cu to the frequency of oscillations. 

Optical cavities with small and high decay rate are required for various applications in op- 
tics [2T| (including spontaneous emission control, low-threshold lasers, low-power switches, and 
adiabatic wavelength conversion) because they generally enhance intrinsically small light-matter 
interactions. Since light is essentially difficult to localize, it is hard to realize small-sized optical 
cavities with strong light confinement. Recently, however, the rapid progress has been made on this 
issue, see ID [13 and references therein. As a result, numerical aspects of emerging optimization 
problems for QN eigenvalues have attracted great current interest among specialists in Applied 
Mathematics HHH1I3- 

An attempt to provide an analytical background for QN eigenvalue optimization problems lead 
the author to the development of the two-parameter perturbation approach [9]. The paper is 
focused on mathematical details of the proofs, and uses techniques of nonlinear functional analy- 
sis and multi-dimensional complex analysis. The present communication is a continuation of the 
paper Q. Its goal is to apply the two-parameter perturbation method to a concrete engineering 
problem and, using this example, to make the main ideas understandable for specialists in applied 
sciences. We consider the problem of optimization of a single QN eigenvalue in one-sided leaky 
optical cavity. While the corresponding eigenvalue problem is derived under several simplifying 
assumptions, it is proved to be very useful and have been extensively studied in Physics and nu- 
merical papers, see e.g. ll20l[T5l [5l. We give rigorous mathematical definitions of locally extremal 
QN eigenvalues and of optimal cavities. In Section [3T21 the existence of such cavities is shown. 
Then the two-parameter perturbation method is introduced in Section @] and applied to derive a 
nonlinear equation for QN modes with extremal properties in Section [61 We use this nonlinear 
equation to reduce the cavity's optimization, which is essentially an infinite-dimensional problem, 
to the four-dimensional problem of finding the zeroes of a specially introduced function W, which 
depends only on one real and one complex variable. The main conclusions are summarized in the 
last section. 

Some of purely mathematical details of the presented method are sketched to make the paper 
more accessible for non-mathematicians. However, a reader interested in rigorous analytical proofs 
will be able to recover them without difficulty using the proofs in [0. 

Note that, in Engineering and Physics literature, the rate of energy loss is often characterized 
by other parameters, like Q-factor or lifetime of a QN mode, which in many cases can be approx- 
imately expressed in terms of Imw and Re a; lfl4ll . Keeping this in mind, it is easy to rewrite the 
results of the paper in terms of these, to some extend equivalent, parameters. 



2 Physical settings and the problem statement 

Consider an isotropic non-dispersive transparent medium which has the relative magnetic per- 
meability equal to 1 everywhere, but which is electrically inhomogeneous in one direction. The 
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electromagnetic structure of this medium is described by the (dielectric) permittivity e, which is 
> 1 and varies only in the x 3 -direction e = e(x 3 ). When e is piecewise constant, this model is 
called a 1-D photonic crystal. 

Under additional simplifying assumptions that electromagnetic waves pass normally and the 
electric field E is in the x 2 -direction E = {0, E, 0}, the study of electromagnetic wave propagation 
can be reduced to the scalar 1-D wave equation 

d^E(xs,t) = ^$E(xa t t), (2.1) 

where the 22-component E of the electric field depends only on x 3 , and c is the speed of light in 
vacuum. 

To introduce a one-sided open optical cavity, assume that the medium has the above structure 
in the slab < x 3 < I with piecewise continuous e (that is, for a certain partition = sq < 
si < • • • < s„_i < s n = I, s(s) is continuous in the intervals s 5 -_i < s < Sj and has one-sided 
limits at points sj). The outer medium is assumed to be a perfect conductor for x 3 < 0, and a 
homogeneous dielectric with constant permittivity > 1 for x 3 > I. The boundary condition for 
the electric field on the plane x 3 = is E = 0. Since the medium is homogeneous for x 3 > /, 
the waves radiated from the cavity are never reflected back. If all sources of waves are inside the 
cavity < x 3 < I, the waves satisfy the radiation boundary condition d X3 E + d t E = at the 
interface x 3 = I. For a monochromatic field E(x 3 ,t) = E(x 3 )e^ lu>t in the cavity, these settings 
lead to the eigenvalue problem 



E»(x 3 ) = -co 2 ^E(x 3 ) (2.2) 
cr 



equipped with the boundary conditions 

E(0) = 0, iu^-EQ) = E'(l). (2.3) 

c 

Complex eigen-parameters cu = a — i/3 are called quasi-normal (QN) eigenvalues. The real 
part a = Re u is the frequency of the standing wave E(x 3 )e~ lult . Since the energy leaks from the 
cavity, the oscillations are decaying. The minus imaginary part /3 — — Imu; is positive for all QN 
eigenvalues and is called the decay rate. 

Let 1 < Ei < 82- Restricting ourselves to cavities with e in the range 

ei < e(x 3 ) < £ 2 , (2.4) 

consider the following optimization problem: devise, for a given frequency range ct\ < a < ol% a 
cavity that has a QN eigenvalue u = a — i{3 with the minimal possible decay rate f3. 
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3 QN eigenvalues: basic properties and extrema's definitions 
3.1 QN eigenvalues and their multiplicities 

Since the QN eigenvalue problem is one-dimensional, from now on we will write s instead of x 3 . 
Let the length / of the cavity and the permittivity of the outer medium s > I be fixed. Then the 
cavity is completely described by the permittivity e(s) in the interval < s < I and, for brevity's 
sake, we will speak on QN eigenvalues of the cavity e. 

By (p(s) = <p(s, z; e) and ip(s) = 4>(s, z; e) we denote the solutions of E"(s) = —z 2 ^-E(s) 
satisfying 

<p(0, z\ e) = d s ip(0, z; e) = 1, d s (p(0, z; e) = ip(0, z; e) = 0. 
In particular, xj) satisfies the integral equation 

rj}{ 8 ) = s - — / (s - t) e{t) ip(r) dr, < s < I. 



c JO 

The QN eigenvalues of e are exactly the roots of the equation 

F(z) = 0, where = F(z;e) =iz^^>ip(l,z;e) - d s ip(l,z;e). (3.1) 

c 

The function F is analytic in z in the whole complex plane C. A nontrivial solution E(s) to the 
problem (12.21 ), (12.31 ) is called a QNmode corresponding to a QN eigenvalue u. Recall that a solution 
E is called trivial if E(s) = for all (more precisely, almost all) s in the interval < s < I. 

All QN modes of e have the form Cip(s, u; e), where C is an arbitrary non-zero constant. Hence 
the geometric multiplicity of any QN eigenvalue u equals 1 . In the following, the multiplicity of tu 
means its algebraic multiplicity. By definition, the algebraic multiplicity of a QN eigenvalue cu is its 
multiplicity as a zero of F |fTT|. The w-dependent boundary condition iu^^-E(l) = E'{1) makes 
the QN eigenvalue problem non-self-adjoint. In particular, QN eigenvalues may be degenerate, 
i.e., of multiplicity > 2, see e.g. [fT3l [6l [T9). Concerning the contemporary theory of spectral 
problems with an eigen-parameter in boundary conditions, we refer to [|2l and references therein. 

Concerning other properties of QN eigenvalues, it is easy to derive that: (i) the multiplicity of 
each QN eigenvalue is finite, (ii) the set of QN eigenvalues is symmetric with respect to (w.r.t.) 
imaginary axis (together with multiplicities), (iii) QN eigenvalues are isolated, oo is their only 
possible accumulation point. The property (ii) allows one to restrict the study of QN eigenvalues 
to the case Re u > 0. 

For a homogeneous cavity with constant permittivity e > 1, the QN eigenvalues are roots of 
tan(— cuy/e l/c) = i^e/Soo. When e ^ e^, they are non-degenerate and form a uniformly spaced 
sequence 

e + e c 
e-e c 

The case when the space outside the cavity is vacuum (i.e., 1 = < e) is discussed in detail in 
ll2~0l . In the case e = e^, there are no QN eigenvalues since the medium is homogeneous for s > 
(the energy of the initial disturbance localized in the region < s < I escapes this region in finite 
time, see e.g. BU for detailed explanations). 



u n = —i——= In 

2/y^ 



TTC j Tl if £ < Eoo 

ViU+lA ife>e OQ ' " = -.-1.0.1..- (3-2) 
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3.2 Local and global extrema for admissible frequencies 

The family Adm of admissible cavities is defined by condition (12.41) . Mathematically, it is con- 
venient to include in Adm all L°° functions satisfying (12.4b . An admissible direction h for an 
admissible cavity e is a function on the interval < s < 1 such that the cavity e + h is admissible 
(we call such h a direction since, for all constants C in the range < C < 1, the cavity e + Ch 
is admissible). We say that cu is an admissible QN eigenvalue if it is a QN eigenvalue of some 
admissible cavity, and say that a is an admissible frequency if a = Re u> for some admissible QN 
eigenvalue cu. 

To get some understanding of admissible frequencies, consider homogeneous cavities. Chang- 
ing constant e in the range E\ < e < e 2 and taking io n given by (13.21 ), we see that the frequen- 
cies Reo; n are admissible. They form a sequence of intervals. For n large enough, the inter- 
val produced by Retu n overlaps with the next one produced by Reu n+ i. Hence frequencies a 
with \a\ large enough are admissible. For instance, let us look closer at the case < e\. 
When constant e runs over the interval E\ < e < e 2 , the real part of lj!hk runs over the interval 

(n + l/2)ej 1/2 < a < [n + 1/2)^ 1/2 . When n > j- j^tjwzi - |. there 

is no gap between the 

n-th and (n + l)-st intervals. Hence a frequency a is admissible whenever 



(here \s] is the smallest integer that is greater or equal to s). 

If an admissible cavity e has a QN eigenvalue u, we say that {cu, e} is an admissible pair. The 
decay rate of an admissible pair is defined by 



While e does not participate in the right hand side, it participates in the domain of definition of the 
functional Dr. This is essential for the next definition. 

We say that an admissible pair {a;, e} is a local minimizer of Dr for a frequency a if, for all 
sufficiently small admissible directions 5e that keep the frequency of a perturbed QN eigenvalue 
lo + 5lo of the cavity e+5e equal to a, the decay rate of {co + 5lo, e + Se}is > Dr(w, e). The rigorous 
form of this condition is: there exist 5 > such that Dr(w, e) > Dr(cj, e) for any admissible pair 
{S3, e} satisfying 



Local maximizers of Dr for a particular frequency are defined in a similar way. 

We say that (3 m i Q (ct) is a minimal decay rate and e is an optimal cavity for a frequency a, if the 
pair {a — i/3 min (a), 5} is admissible and /3 m - m (a) < Dr(o3, e) for any admissible pair {u3, e} with 
Re ui = a. In other words, such a pair {a — i/3 m i n (a;), e} is a global minimizer (over the admissible 
family Adm) of the decay rate for a frequency a . We will also say that the corresponding QN 
eigenvalue a — i/3 m in(ci) is optimal for the frequency a. 




Dr(w,e) 



Imcu. 



e(s) — e(s)| < 5 for almost all < s < I, 
\uj — uj\ < 5, and Re cu = a. 
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Lemma 1. For each admissible frequency, there exists an optimal cavity. 

The existence of optimal cavities follows from the fact that if admissible QN eigenvalues u n 
approach a complex number to, then to is an admissible QN eigenvalue. To explain this fact, 
consider admissible cavities e' n ' that produce the QN eigenvalues ui n . Using the Gronwall-Bellman 
inequality and assumption ( 12.41 ), it is easy to show that the sequence of related ^-solutions ip n {s) = 
ip(s, u n ; £") is uniformly bounded in s and n together with the sequences of their s-derivatives ip' n 
and i\)" n . Hence some subsequence ip nj converges uniformly to a continuous function ip(s) (due to 
the compact embedding of V^ 2,00 ^ 1] into C[0, 1]). Then applying the Banach-Alaoglu theorem 
for L°° (0, 1) to the sequence e^ Tlj \ one can pass to a limit in the integral reformulation of (I2.2I) - (I2.3I) 
(see BU, and also lfT2l for compactness arguments applied to a self-adjoint eigenvalue optimization 
problem). 



4 Two-parameter perturbations 

We study properties of local extremizers (i.e., local minimizers and maximizers) via the pertur- 
bation theory for QN eigenvalues using the following arguments. Assume that {co,e} is a local 
minimizer. If we consider all small admissible directions Se that keep Re(co + Stu) equal to a, 
then corresponding corrections 5(3 = — Im Su of the decay rate must be positive. The key point is 
to show that there are enough various 5e with the above property to impose strong restrictions on 
extremal e. 

For non-degenerate QN eigenvalues, the first order perturbation correction can be calculated 
by the standard sensitivity analysis. However our method requires a more delicate approach due 
to following reasons: (i) QN eigenvalues may be degenerate, (ii) we are interested only in the 
perturbed QN eigenvalues that stay on the line Recu = a. Indeed, even if one is able to find 
the first order correction that moves a QN eigenvalue along the line Re a; = a using singular 
perturbation theory (see e.g. |[T6ll and references therein), there are no any guarantee that higher 
order corrections does not contribute to the real part of oj. 

To resolve these difficulties we study two-parameter perturbations of QN eigenvalues consider- 
ing them as z-roots of the equation F(z, e + Se) = 0. This requires computation of the (complex) 
partial derivative d z F(z; e) and the directional derivatives d^F(z; e + (h) for admissible pairs 
{a;, e}. Here and below £ is a complex number, and h is an arbitrary (measurable and essentially 
bounded) function on < s < I, which is called a direction. 

The derivative d^F(z; e + (h) exists for an arbitrary direction. More precisely, the Maclaurin 
series for the solution ip (see e.g. 0) can be used to get the Maclaurin series for F(z; e) and, in 
turn, to show that for arbitrary n directions hi(s), . . . , h n (s), the function F(u + z\ e + Q\hi + 
^2^2 + • • • + Cnh n ) is an analytic function of n+1 complex variables z, 0, . . . , ( n (moreover, the 
map {z, s} h-> F(z; e) is analytic on C x L°°(0, 1)). 

Differentiating the integral equation for tp and solving the produced differential equations, it is 
possible to express d z F(z; e) and d^F(z; e + (h) via the fundamental solutions if and if). When 
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z = cu is a QN eigenvalue of e, the obtained expressions can be simplified to 

d ( F(u;e + (h) = — -f ^ 2 (s,cu;e)h(s)ds, (4.1) 

d z F{u-e) = — r / ij 2 (s,u;e)e(s)ds + i^iP(l,u;e). (4.2) 

c 2 ip{l,cu;e) J c 

The last equality yields the following lemma. 

Lemma 2. A QN eigenvalue cu of e is degenerate exactly when 

' C\/E 
i) 2 (s,u;e)e(s)ds + i y °° 4> 2 (l,cu;e) = 0. 

leu 



Note that 



VjJ- 



-<p(l,u;e) - tff(l,u\e) 



7^0. 



Let us consider the first-order correction for cu in the case when cu is a non-degenerate QN 
eigenvalue of e. Small perturbations of e lead to small perturbations of cu and the perturbed QN 
eigenvalues remains non-degenerate. Taking 5e = Qh with a fixed perturbation direction h = h(s), 
it is easy to show that the perturbed QN eigenvalue cu + 5cu is an analytic function fi(C) of ( in a 
vicinity of Co = 0. The approximation ofQ(Q to the first-order term is given by 

, , w fn4> 2 (s)h(s)ds , „ x 

n(C) = w-c — i V- + °C ■ (4.3) 

2/>2 (sMs)ds + i ^ 2(/) 

This formula, as well as higher order corrections in the non-degenerate case, are known [15]. 

Consider the case when cu is degenerate. Perturbations of a multiple eigenvalue may lead to its 
splitting in several eigenvalues. The general splitting picture is quite complicated and is described 
by one or several series (Puiseux series) in fractional powers of ( lfT8lfT0lfT6l . For our needs, it is 
enough to consider only the generic case when the direction h satisfies d^F(cu; e + (h) ^ . This 
condition is equivalent to 

/ ip 2 {s,cu;e)h(s)ds ^ 0. (4.4) 
Jo 

Assume that cu is a QN eigenvalue of e of multiplicity m, and that perturbation's direction h 
satisfies (14.41) . To study the equation F(cu + z\ e + Qh) = 0, we introduce the function Q(z, C) = 
F(cu + z; e + Qh) and write it as a power series in z and £. The equation Q = takes the form 

d m O(0 0) 

z m z v \ ' } + c9 c g(o, o) + o(z m ) + o(0 = o 

ml 

for z and ( going to 0. Note that d™Q(Q, 0) ^ since cu has multiplicity m. Considering z as a 
function of £, we see that |z(C)| < C|C|^ m for C small enough with some constant C. So o(z m ) is 
also o(C). Hence, 

m!g c g(o,o) 
lu ; flfg(o,o) + lu ' 
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According to ( 14.11) and the assumption on h, 

,2 r l 



d c Q(0,0) = —£ r / ^ 2 ( S , W ;e)/i( S )d^0. 

c 2 i;{l,uj;e) J 

Taking the m-th root, we get a formal first-order correction formula for a perturbed QN eigenvalue 

n(C) 

n(C)=w+(tCi(u,e,m) ^(s,u;e)h(s)ds) +o(|C| 1/m ), (4.5) 



where 



m\uj 2 



C 1 (w,e,m) = — — rTT-rr, c- (4.6) 

c 2 ip(l, u; e)ofF(<jj, e) 

The m-th derivative d™F{u, e) can also be expressed through the fundamental solutions if and ip, 
but we are not going to use such a refinement in this paper. Note that (14.51) turns into (14.31) for 
m — 1. 

Now consider two-parameter perturbations of u, assuming that the directions h\ and h 2 are 
such that d^F(z; e + Ci^i) and d^ 2 F(z; e + C2^2) are nonzero and 

0< kTgd^F{z-e + C2h 2 ) - A-Tgd Cl F{ Z] e + Cihi) < tt (4.7) 

for a suitable branch Arg of the arg-function (here and below arg z is the complex argument 
or phase of a nonzero complex number z). Condition (14.71) means that d^F(z;e + Ci^i) and 
d{ 2 F(z; e + (2I12), perceived as vectors must not be oppositely directed and, additionally, must be 
ordered such that (14.71) holds. Considering a perturbed eigenvalue u + z of e + Ci^i + C2^2 as a 
function of £i and £ 2 > we write the equation 

Q(z,G,C 2 ) = 

for the correction z = z(Ci, C2) defining Q as 

Q{z, Ci, C2) = F(u + z, e + Cihi + <M). 

Denote C = {&, C2}, |C| = (Ci + C 2 2 ) 1/2 > and = {0, 0, 0}. When |C| is small, we have 

Z -" Cl d?Q0 ~ Ca d?Q(0) + ° {Z ) + °(C0 + o(C 2 )- (4-8) 

Assume additionally that 

(1 and (2 are nonnegative and of the same order. (4.9) 

That is, (1 and ( 2 go to zero such that their ratios d / £ 2 and £ 2 / Ci are bounded. Under this condition, 
( Ob and (S3 yield that (2)! < c|C| 1/m for small |C|. Hence, 

/ m!%Q(0) m!%g(0) \ 1/m , 
*(Ci, C2) - l^-Ci dTQ{Q) - C2 dTQ{Q) J +o(\(\ ). (4.10) 
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Appending the values of partial derivatives provided by (14.11 ), one gets the first-order approxima- 
tion for perturbed QN eigenvalues 

/ a \ 1 / m 

n(C 1 ,C2)=u}+\C 1 {u,e i m) J ^ 2 (s,cu;e)[C 1 h 1 (s) + C2h 2 (s)}ds\ +o(|C| 1/m ) (4.11) 

with Ci defined by (|4!6t . 

For sufficiently small positive Ci> C2 satisfying ( 14.91) , the m branches of m-th root produce m 
different perturbed QN eigenvalues, each of which is non-degenerate. The proof of this fact can 
be given along the lines suggested in the paper flU, where a slightly different eigenvalue problem 
was considered. Condition (14.71 ) ensures that the coefficient under the m-th root is nonzero. It is 
important for the sequel that this coefficient is a linear functional of perturbation directions hi and 
h 2 . 

Formula (14.1 II) will be applied to the study of the QN eigenvalue optimization via the lemma 
in the following section. 



5 Admissible directional derivatives and arguments of eigen- 
value corrections 

Let us fix an admissible pair {u, e}. The set of admissible directional derivatives of F at {u, e} 
is defined as the set formed by complex numbers d^F(ui; e + (h) with h running through the set 
of admissible directions (i.e., we take all h such that e + h is an admissible cavity). The set of 
admissible directional derivatives contains and is convex (that is, if £1 and £2 are admissible 
directional derivatives, then 7^1 + (1 — 7)^2 is so for all < 7 < 1). We say that the set 
of admissible directional derivatives contains a neighborhood of zero if it contains all complex 
numbers £ with modulus |£| less than a certain positive number. Let Arg be the branch of the 
arg-function that takes values in the interval — ir < Arg z < n. 

Lemma 3. Assume that the set of admissible directional derivatives of F at an admissible pair 
{u, e} contains a neighborhood of zero. Then for any 9 in the range —it < 6 < ix, there exists a 
sequence of admissible pairs {co n , e^} such that 

(i) Arg (oj n — u) = 6 for all n = 1,2,..., 

( ii) co n tends to to and tends to e ( the latter means that the essential suprema of\e(s) — e^ n \s)\ 
tends to zero as n goes to 00 ). 

In particular, the lemma yields that if{cu, e} is a local extremizer ofDifor the frequency Re u, 
than the set of admissible directional derivatives of F at {u, e} cannot contain a neighborhood of 
zero. 

To derive the lemma, note that its assumptions allows us to chose admissible directions hi and 
h 2 and the branch of the m-th root in (14.1 II ) such that 
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(Al) assumption (14.71) is valid for hi and h 2 , 

(A2) for positive £, the corresponding first order corrections to u, 

i v i/m / r l \ 

(Ci(uj,e,m) / ?/; 2 (s, 5)/ii(s)ds J and l£Ci(u;,£, m) / ip 2 (s, u; e)h 2 (s)ds J , 

have the arguments slightly smaller and, respectively, slightly greater than 6, 

(A3) the chosen branch of the m-th root is continuous in the sector of complex plane covered by 
numbers ^ 

Ci(cu,e,m) / ^ 2 (s,cu;e)[Ci/ii(s) + ( 2 h 2 (s)]ds 
Jo 

while d and £ 2 run through positive numbers. 



This allows us to apply formula (14. lib , which yields that the argument Arg (fi(Ci, C2) _ of 
the QN eigenvalue's perturbation changes continuously for small positive d and £ 2 . This and 
(A2) imply that for small enough e there exist positive £i and ( 2 such that Ci + C2 = e and 
Arg (Jl(Ci> C2) — equals (up to a multiple of 2ix). Finally, we take such £1, ( 2 for a sequence e n 
going to 0. Note that the convexity of the families of admissible cavities and admissible directions 
is essentially used here. 

6 Nonlinear eigenvalue problem for local extrema 

In this section we derive a nonlinear eigenvalue problem for QN modes corresponding to local 
extremizers. 

Lemma 4. Let e be an admissible cavity and \mz 2 < 0. Then ip(s,z;e) 7^ in the interval 
< s < I, and the function £(s) = argip 2 (s, z; e) is strictly increasing (the values of sag are 
chosen such that £ is continuous for < s < I). 

Let us show this. Assume first that ip(si) = 0, where ip(s) = ip(s, z\ e). Then 

z 2 r \ip(s)\ 2 e(s)ds = -c 2 f * ip"(s)i^s)ds = -c 2 H \f(s)\ 2 ds. 
Jo Jo Jo 

The number in the right-hand side is real, but by our assumptions z 2 is not real. Hence s± = 
(note that ip'(Q) = 1 and so tp is non-zero for small positive s). Thus, tp(s) 7^ for s > 0. Note 
that 

^) = 2I m (ln^( S )y = 2Im^^. 

|^(s)| 2 

This shows that, on the one hand, the function \ip\ 2 ^' tends to as s goes to + , and on the other 
hand, 

(\^(s)\ 2 £(s))' = 2ImV>"(s)#0 = -2c~ 2 e(s)\ij(s)\ 2 lmz 2 > 0. 
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Hence \ip\ 2 i' is positive for < s < I, and so is 

Denote by Xa,b the function that equals 1 in the interval a < s < b and equals outside this 
interval. For a finite number of complex numbers Z\, z 2 , z n , the convex hull of these numbers, 
by definition, consists of numbers C\Z\ + c 2 z 2 + . . . c n z n , where c\ + C2 + ■ ■ ■ + c n = 1 and all Cj 
are nonnegative. Clearly, a set of complex numbers is convex if and only if it contains a convex 
hull of every its finite subset. 

Assume that e is a piecewise continuous admissible cavity, and that u> with Re u > is a QN 
eigenvalue of e. Since changing of e at a finite number of points does not influence QN eigenvalues, 
this is not an essential restriction to assume from now on that e takes its limit from the left or its 
limit from the right at the points of discontinuity and the endpoints s = and s = I. 

Further, suppose that e\ < e(s) < e 2 in a certain interval si < s < s 2 . Then 

E\ + 8\ < e(s) < e 2 — 5\ in some narrower interval s 3 < s < s 4 (6.1) 

for small enough Si > 0. Hence, the two directions §iXa,b and (— l)Six a ,b are admissible for every 
a and b such that s 3 < a < b < s 4 . It follows from Lemma H] that there exist ai,a 2 , b±, b 2 with the 
properties: 

(i) s 3 < ai < bi < a 2 < b 2 < s 4 , 

(ii) the directional derivatives d(F(z;e + CXoi,6 1 )> d^F{z;e + CXa 2 ,b 2 ) are non-zero and the 
difference of their arguments is not a multiple of n. 

With such a choice of a\ 2 and b X 2 , the convex hull of the four admissible directional derivatives 
d(F{z\£ ± C^iX<n,6i)> dc t F(z]S ± (5iXa 2 ,b 2 ) contains a neighborhood of zero. Since the set of 
admissible directional derivatives is convex, it also contains a neighborhood of zero. By Lemma[3] 
the pair {u, e} is not a local extremizer. 

Thus, we have shown that if an admissible pair {cu, e} (with piece-wise continuous e) is a local 
extremizer, then e takes only the extreme possible values E\ and e 2 . 

Since the set of admissible directional derivatives at {tu, e} is convex, it either contains a neigh- 
borhood of zero or, otherwise, 

the set of nonzero admissible directional derivatives is contained in 

a certain complex half-plane 9 < arg z < 9 + ix. (6.2) 

Assume again that {u, e} is a local extremizer and e is piecewise continuous. Then Lemma[3] 
yields that the case (16.21 ) takes place. We have shown before that e takes only values E\ and e 2 . 
Now, repeating the use of Lemma [3j we want to find restrictions on the points s where such an 
extreme e switches between e\ and e 2 . To this purpose, consider directional derivatives produced 
by directions ±5±x a ,b with intervals a < s < b of small lengths. For any point so, there exists an 
interval (of positive length) s 4 < s < s 2 containing s and such that e is constant on it. To be 
specific, assume e(s ) = S\. Then for any subinterval a < s < b of the interval s x < s < s 2 , 
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the direction 5\Xa,b is admissible whenever < 8\ < e 2 — E\. The corresponding admissible 
directional derivative d^F(u, e + C^iXa,b) equals 



u 2 /■» 



c^(f,w;e) J a 

Taking a and b close enough to s , it is possible to make the difference of arguments 

&Ygd ( F(uj,e + (5!Xa,b) - arg — — r^ 2 (s , u; e] 



c 2 ^(/, uj; e) 

arbitrary small (after possible correction on a multiple of 2ir). Combining this with (16.21 ), one can 
see that 

#o < arg ( -r—TT- r^ 2 (s , w; e) J < O + tt. 

Similarly, if e(so) = £2, then there exist a and b arbitrary close to so such that the direction 
(— l)SiXa,b is admissible whenever < 61 <e% — E\. This yields 

O - tt < arg f — — ^ r^ 2 (-s , w; e)^ < ^o- 

Summarizing, we see that, for a local extremizer with piece-wise continuous e, there 

exists Q\ such that e takes the value S\ in the intervals where B\ < &rgip 2 (s, oj\ e) < 61 + n, and 
takes the value £2 in the intervals where 9i—n< arg i(j 2 (s, u; e) < 61 (here, the double inequalities 
with multi-valued arg function are assumed to be valid if they are valid for one of its values). The 
values of e at the finite set of points where ^ 2 (s) cross the lines arg z — Q\ and arg z = 61 + tt are 
not important since these values do not influence QN eigenvalues of e. The angle 61 is related to 
9 of dO) through 

^ = 00 -arg— -. (6.3) 

Similar arguments also works for the case Re to < (for instance, trough the use of the sym- 
metry of QN eigenvalues w.r.t. iW) and for the case Re cu = (for this case some modifications are 
needed, but considerations are simpler since the solution ip is real-valued). Finally, note that the 
assumption of piecewise continuity of e can be dropped (that is, e can be assumed only to be an 
admissible L°° -function). This can be done if instead of small intervals a < s < b, one considers 
sets of small measure and small diameter. Corresponding modifications require a number of purely 
technical details, which a concerned reader can recover without difficulties using 

It is convenient to write the obtained result in terms of a nonlinear eigenvalue problem. Denote 
by {Imz > 0} the (open) upper complex half-plane, which, by definition, consists of complex 
numbers z with positive imaginary part. Similarly the half -plane {Im z < 0} consists of numbers 
z such that lm.z < 0. Let Xi m *> (C) = 1 when Im^ > 0, and Xi mz >o(C) = when Im^ < 0. 
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Theorem 5. Assume that {oj,e} is a local extremizer (i.e., minimizer or maximizer) of the de- 
cay rate for an admissible frequency a = Re u. Then the boundary value problem ( 12.31 ) for the 
nonlinear equation 

E" = ~E [e, + (e a - e 1 )x Img>0 (^)] (6.4) 
has a nontrivial solution E(s) satisfying 

e(s) — E\ + (e 2 — £ i)Xi mz>0 (E { s )) (almost everywhere) on < s < I. (6.5) 

Let us show this for the case Recu > 0. Multiplying the solution ip to a constant e^"" -01 " 2 
with 9\ from (16.31) . one gets another solution E(s) = e 1 ^ 6l '^ 2 il>(s, u>; e) to the linear equation 
E" = —^sE. This solution additionally satisfies (12.31 ) and has the property that E 2 takes values 
in {Imz < 0} when e(s) equals e l5 and takes values in {Imz > 0} when e(s) = e 2 (after a 
possible correction of e at a finite number of points). So e and E are additionally connected by 
(16.51) and, therefore, E is a solution to (16.41) . Modifications of this proof for the cases Re a; < 
and Re to = are similar to that of BH. 

Note that if y(s) is a solution to the equation 

E" = ~E [5i + (s 2 - ei)x Imz>0 (^ 2 )] , (6-6) 

then Cy(s) is so for each positive constant C. This implies that any nontrivial solution y to ( 16.61 ) 
satisfying y(0) = may be written in the form y(s) = C^/(s; z, 9) with a positive constant C and 
a function \P(s; 2, 6 1 ) defined by 

a s 2 v]/ = -^M/ [ £l + ( £2 - £l ) Xim ^ >o (^)] ; ^( O;z ,0)=O, 9 a W(O;z,0) = e w , 

where 2; is a complex number and 9 is a number in the interval — ir < 9 < n. It is easy to see that, 
for the initial data E(0), E'(0) such that E(0) = and E'(0) ^ 0, equation (|6T6b has a unique 
solution, and so \& is well-defined. 

For —ix<9<ix and a complex number 2, define the function 

W(z, 9) = iz^^^!(l] z, 9) - d s V(l; z, 9). 
c 

One can see that the nonlinear eigenvalue problem (16.41) . (12.31) has a nontrivial solution exactly for 
those numbers cu that satisfy 

W(u, 9) = for at least one 9. (6.7) 

Let (16.71 ) be valid with a certain 9. Then y(s) = ^(s; w, 0) is a solution to the original linear 
eigenvalue problem (12.21 ). (12.31 ) with admissible s(s) = S\ + (e 2 — 5i)x Imz>0 (z/ 2 ) ? and w is a QN 
eigenvalue of e. Thus, w is an admissible QN eigenvalue. This, Theorem [5J and the definitions of 
optimal cavities and QN eigenvalues (see Sect J3.2l) yield our main result. 
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Theorem 6. The minimal decay rate /3 m i n («) for an admissible frequency a equals to the minimal 
real number (3 having the property that W(a — 9) = for at least one 9 in the range —tt < 

9 < 7T. 

Moreover, let uj = a — i/3 m in(«) be the corresponding optimal QN eigenvalue. Then for any 9 
such that W(uj, 9) = 0, the cavity defined by 

e(s) = ei + (e 2 - £i)x Imz>0 {^(s; u, 9)) (6.8) 

is optimal for the frequency a. And vice versa, for each cavity e optimal for a , there exists 9 such 
that d<5.g|) holds. 

7 Conclusions 

The paper is concerned with QN eigenvalues of 1-D leaking optical cavities. The cavities are 
described by dielectric permittivity function e which depends on one variable 23 and is assumed 
to take values in a fixed range £1 < e(x-s) < £2- Cavities satisfying this restriction are called 
admissible. 

We study analytically cavities that produce QN eigenvalues cu with a given frequency Re u = a 
and locally extremal (locally maximal or minimal) for this frequency decay rate (3 = -Imw. We 
show that cavities with such locally extremal properties are 1-D photonic crystals consisting of 
alternating layers of two materials with two extreme allowed dielectric permittivities e\ and £2- 
This explains effects observed in numerical experiments for very kindred (but slightly different) 
optimization problems (81 13. To find thicknesses of the layers in extremal cavities, we derive 
a nonlinear eigenvalue problem for their QN modes. It occurs that x 3 coordinates of interface 
planes between the layers are tied to rotation in the complex plane (i.e., to the arg-function) of the 
corresponding extremal QN mode. 

For each admissible frequency a there exist a cavity that creates a QN eigenvalue with minimal 
possible decay rate. In the paper, this QN eigenvalue is called optimal for the frequency a. We 
show that such optimal QN eigenvalues can be easily found via zeroes of a function W, which 
is constructed by solutions of the above mentioned nonlinear equation and depends only on two 
variables. This effectively excludes the unknown optimal dielectric permittivity e from the process 
of calculation of an optimal QN eigenvalue. After calculation of an optimal QN eigenvalue and 
corresponding QN modes, optimal dielectric permittivity functions can be easily recovered from 
their connection with optimal QN modes. 
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